Guided Wave Propagation Study on Laminated Composites 
by Frequency-Wavenumber Technique 

Zhenhua Tian ', Lingyu Yu '*, Cara A.C. Leckey 2 
'Mechanical Engineering Department, University of South Carolina, Columbia, SC 
Nondestructive Evaluation Sciences Branch, NASA Langley Research Center, Hampton, VA 


ABSTRACT 

Toward the goal of delamination detection and quantification in laminated composites, this paper examines guided 
wave propagation and wave interaction with delamination damage in laminated carbon fiber reinforced polymer (CFRP) 
composites using frequency-wavenumber (f-k) analysis. Three-dimensional elastodynamic finite integration technique 
(EFIT) is used to acquire simulated time-space wavefields for a CFRP composite. The time-space wavefields show 
trapped waves in the delamination region. To unveil the wave propagation physics, the time-space wavefields are further 
analyzed by using two-dimensional (2D) Fourier transforms (FT). In the analysis results, new f-k components are 
observed when the incident guided waves interact with the delamination damage. These new f-k components in the 
simulations are experimentally verified through data obtained from scanning laser Doppler vibrometer (SLDV) tests. By 
filtering the new f-k components, delamination damage is detected and quantified. 
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1. INTRODUCTION 

Current composite “safety margin” design comes with a weight penalty which could be mitigated by the 
development of accurate, efficient, and reliable structural health monitoring (SHM) methods to locate and quantify the 
unique damage types that can occur in composite materials. Unexpected damage can occur in aerospace composites due 
to impact events or due to stressing of the material during off-nominal loading events. Matrix cracking, fiber breakage, 
fiber-matrix debonding, and delamination are the most common types of damage in fiber composites [1]. In particular, 
laminated composites are especially susceptible to delamination damage due to their weak transverse tensile and 
interlaminar shear strengths (as compared to their in-plane properties)[l]. Delamination damage can occur internally 
such that it is not readily visible on the composite surface [2]. Hidden delamination damage caused by the impacts must 
be detected and evaluated before it becomes critical. The problem of how best to detect such material changes in 
laminated composite components is still being solved, and a solution to this problem will be imperative for safe and 
functional, optimally-designed, next-generation composite aircraft. 

Guided ultrasonic waves have proved effective for SHM due to their ability to inspect large areas while 
maintaining sensitivity to small defects in the structure; and also due to the capability to implement guided wave 
techniques using permanently installed, lightweight low-profile piezoelectric wafer sensors [3-9]. Advances in guided 
wave based SHM technologies over the last decade have demonstrated the feasibility of detecting and locating damage 
in composite structural components [10-19]. Due to the complex wave interaction that occurs when hidden delamination 
damage is present, tremendous efforts have been put into delamination detection methodologies. Guo and Cawley 
studied the interaction of the S 0 Lamb mode with delaminations [20], showing that the reflection amplitude of the S 0 
mode from a delamination is strongly dependent on the delamination depth. Petculescu et al. used a group delay 
measurement technique for the detection and sizing of delaminations in unidirectional and cross-ply composites [21]. 
They found that the delamination affects the traveling time of ultrasonic wave packet. Staszewski et al. used 3-D laser 
vibrometry to locate delamination damage and estimate its severity in a composite plate, demonstrating that the 
delamination can be revealed by the amplitude profiles of Lamb waves [2]. Purekar and Pines demonstrated the 
application of piezoelectric phased sensor arrays to detect the presence of delamination damage in a composite laminate 
[22]. Michaels et al. used guided wavefield images and frequency- wavenumber domain analysis to study the wave 
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interactions with structural discontinuities in composites [23]. Yeum et al. presented a delamination detection technique 
using a dual piezoelectric transducer network [24]. The result for the case they studied showed that the anti-symmetric 
A 0 mode slowed down when it passed through a delamination region, while the speed of the symmetric S 0 mode was 
little affected by the delamination. 

All of this prior work demonstrated that delamination damage affects various characteristics of the interrogation 
waves (such as amplitude and travelling time). In addition, the damage can cause wave scattering, mode conversion, and 
multiple reflections in the delamination region. To understand the mechanisms of wave propagation and interactions in 
the delamination region, other researchers studied the mechanisms through simulation and experiment. Hayashi and 
Kawashima reported the use of numerical simulations to study guided wave propagations and interactions in delaminated 
regions [25]. Multiple reflections in the delamination area were confirmed. Hence, a considerable amount of ultrasonic 
energy was “trapped” above and below delaminated regions until eventually being dissipated. Ramadas et al. studied the 
interaction of the anti- symmetric A 0 mode with symmetric delaminations [26]. They also identified multiple reflections 
in the delamination region. Sohn et al. have studied the trapped phenomenon of guided wave propagation in delaminated 
composites through the use of wavefield images obtained from a SLDV[27, 28]. It was observed that upon entering the 
delaminated region, the incident waves would split above and below the delamination and then propagate independently 
through upper and lower lamina (single delamination scenario being assumed). After that a significant portion of these 
waves were reflected back from the far edge (exit region). The reflected waves then pass back above (or below) the 
delaminated area and experience reflections at the original entrance. The process was repeated until the energy was 
dissipated, as found in the simulation work reported by Hayashi and Kawashima [25]. Glushkov et al. also observed 
similar wave energy trapping and localization in a strip delamination by using laser vibrometry to record the wavefield 
[29]. Despite these advances, guided wave techniques for composites are still insufficient in regards to a quantitative 
understanding of the mechanisms of wave interactions with delaminations in composites. 

In this paper, we examine guided wave propagation and interaction in the delamination region in the carbon fiber 
reinforced polymer (CFRP) composites by using f-k analysis for the purpose of delamination detection and 
quantification. Elastodynamic finite integration technique (EFIT) simulations and SLDV experiment are both used to 
obtain the time-space wavefields regarding the wave interactions with the delamination damage. The f-k analysis is 
implemented to expand the level of details such as f-k components and wave mode contents of the wavefields. The f-k 
filtering technique is used to detect and quantify the delamination. This paper is organized as follows. Section 2 presents 
the guided wave simulation technique that was implemented to acquire in-plane and out-of-plane wavefields, and 
compares the wavefields between the pristine and delaminated plates. Section 3 describes the f-k analysis technique 
which was used to analyze the wavefields, and compares the analysis results between the pristine and delaminated plates. 
Section 4 presents experimental tests to verify the wavefields and f-k analysis results of the EFIT data. Section 5 presents 
the f-k filtering which is used for delamination detection and quantification. A summary of findings is provided in 
Section 6. 


2. EFIT WAVE MOTION SIMULATION 

Computational wave simulation has proven to be an effective tool for clarifying complex ultrasonic wave 
phenomenon without labor-consuming and costly experimental exploration. Previous work by the authors has 
demonstrated how 3D simulated guided wave data can aid in the study of guided wave interaction with damage in 
isotropic and anisotropic materials [19, 30, 31]. In this work, 3D EFIT was implemented to acquire guided wave 
propagation data for laminated composites, including cases of wave interaction with the delamination damage. For 
composite materials, EFIT can incorporate ply-level material property changes; thus allowing for specific CFRP layups 
to be accurately simulated. The code is custom parallelized to run efficiently on cluster and multicore computing 
resources. Additional details about the anisotropic EFIT simulation code, including results of code validation studies for 
pristine composite plates, can be found in prior work [19]. The complete EFIT equations can be found in reference [32]. 
Note that for the EFIT model implemented for this paper, 3D refers to the three Cartesian spatial dimensions, and that 
the wave propagation data generated by the simulation are 4D including time. 
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Figure 1. EFIT simulation configuration: (a) the simulation layout; (b) the side view of the [0 2 /90 2 ] s layup. 




Figure 2. Dispersion curves of the 8-ply IM7/8552 composite laminate with [0 2 /90 2 ] s layup: (a) f-k curves for 90° 
direction; (b) group velocity curves for 90° direction. 


Table 1 Material properties 


p (kg/m 3 ) 

E x (GPa) 

E 2 (GPa) 

E 3 (GPa) 

Gi 2 (GPa) Gi 3 (GPa) G 23 (GPa) Vi 2 Vi 3 v 23 

Thickness (mm) 

1570 

171.4 

9.08 

9.08 

5.29 5.29 2.80 0.32 0.32 0.5 

0.1068 


2.1 EFIT simulation setup 

In this work wave propagation in an 8-ply IM7/8552 composite laminate plate with layup [0 2 /90 2 ] s (Figure lb) is 
studied. The total thickness of the simulated laminate is 0.854 mm. The material properties for a single ply are given in 
Table 1, where E x is in the fiber direction (and where v is Poisson’s ratio). With the composite layup and material 
properties, the dispersion curves along 90° (corresponding to the direction along the dotted line in Figure la) were 
calculated using DISPERSE software [33]. The calculated dispersion curves are plotted in Figure 2. The dispersion 
curves show that below the cutoff frequency 1 MHz there are two guided wave modes existing in the plate, the 
fundamental A 0 mode (with lower velocity) and S 0 mode (with higher velocity). 

Both the pristine and damaged plates were modeled by using the EFIT. To make a simple case for preliminary 
study, the damage is a 20 mm x 20 mm delamination between the uppermost 0/90 plies as illustrated in Figure lb. Prior 
work by Leckey et al. has demonstrated that delamination geometries can be quite complex (depending on the composite 
layup and damage severity), and the simple damage case studied in this work is considered to be a starting point for 
developing a quantitative understanding of in-plane and out-of-plane guided wave interaction with delaminations in 
composites [19, 34]. 

In order to accurately capture the small wavelengths of wave modes that may exist in the layers (above or below 
the delamination as illustrated in Figure lb) which will have much smaller thickness, a small spatial step size was 
enforced for the EFIT simulation. The spatial step size was set to Ax=20 pm for both the pristine and delaminated cases. 


It is important to note that EFIT requires a spatial step size no larger than Ax<^ in /10, and for cases where delaminations 
(or a change in thicknesses) are present, a ^ min based on the full plate thickness may not suffice. The delamination was 
incorporated by implementing stress free boundary conditions at the delamination location. Excitation used in the model 
is a 300 kHz 3 -count Hanning window smoothed sine wave. The excitation was introduced to the simulation space as a 7 
mm diameter ring-type excitation on the top laminate surface in order to approximate the behavior of a 7 mm bonded 
piezoelectric wafer [30]. Velocity data as a function of time (£) and space (x, y) from the top surface of the composite 
plates, including in-plane v* and v y components as well as out-of-plane v z component, are generated from the EFIT 
simulation. 



(a) Time (ps) 




Figure 3. The time-space wavefields of EFIT in-plane v v components: (a) for the pristine plate; (b) for the plate with a 
delamination. 



Figure 4. The time-space wavefields of EFIT out-of-plane v z components: (a) for the pristine plate; (b) for the plate 
with a delamination. The out-of-plane component of A 0 is prevailing compared to that of S 0 . 

2.2 Simulation results 

To study the wave interaction with the delamination, the EFIT simulation results along the y axis from the source 
through the delaminated region (dotted line in Figure la) are used. The selected simulation results are two dimensional 
(2D) in terms of argument y (propagation distance) and argument t (time), denoted as time-space wavefield v(t , y). The 
v(t, y) notation refers to velocity wavefield as either in-plane, v x (t, y) and v y (t, y ), or out-of-plane, v z (t, y ), components. 
The in-plane v* component is zero in this study due to the fact that the v* component is normal to the wave propagation 
direction along the y axis (90° direction). Hence, only the in-plane v y and out-of-plane v z components will be analyzed 
and discussed in this paper. 

Figure 3a and b show the time-space wavefields of the in-plane v y components for the pristine and delaminated 
plates, respectively. A faster S 0 mode and a slower A 0 mode, as well as their interactions with the structural 
discontinuities (plate boundary and delamination for the delaminated cases), can be clearly viewed. The in-plane S 0 
mode shows higher amplitude than that of the A 0 mode. It is seen from the wavefields that in the delamination region, 




waves bounce back and forth between the boundaries of the delamination, distinctively exhibiting the “trapped waves' 
phenomenon. The observed behavior agrees with the wavefield observations presented by Sohn et al [27]. 


Figure 4a and b show the time-space wavefields of the out-of-plane v Z9 components for the pristine and delaminated 
plates, respectively. It is seen the out-of-plane A 0 mode has higher amplitude than the S 0 mode. The “trapped waves” 
phenomenon is also observed in the out-of-plane wave interactions with the boundaries of the delamination, similar to 
that of the in-plane v y observations. 

3. TRAPPED WAVE MODES WITHIN DELAMINATIONS 

In the above section, the time-space wavefield shows wave interactions with the delamination and exhibits the 
trapped waves in the delamination region. However, implicit wave propagation characteristics such as f-k components 
and wave mode contents of the time-space wavefield are not readily seen. It has been demonstrated that wavefield 
analysis representing the f-k content has abundant information regarding the existence of various wave modes and wave 
propagation characteristics [23, 35-37]. Thus, in this section we transform the time-space wavefield to the f-k 
representation to analyze the wave interactions with the delamination. 

3.1 F-k analysis 

The f-k representation can be obtained from the 2D time-space wavefield data by applying the 2D FT, 
mathematically given as [23, 35-38] 


U(f,k) = F 2D [w(7,x)] = J J u(t,x)e j ^ 2nft ^dtdx (1) 

where U(f. ] k) is the resulting f-k representation or “spectrum” in terms of the frequency variable / and the wavenumber 
variable k. The U(f. ] k) can be interpreted as an alternative representation of the time-space wavefield u(t , x). The f-k 
spectrum not only shows the f-k components of the wavefield but also provides an easy means to identify the wave mode 
compositions. 

3.2 In-plane component v y analysis 

The f-k spectrum of the EFIT in-plane v y wavefield for the pristine plate is first calculated and plotted in Figure 55 a. 
Theoretical dispersion curves for waves in the 90° direction in the [0 2 /90 2 ] s laminate are plotted in the same figure as 
dotted (A 0 ) and solid (S 0 ) lines. The pristine plate f-k spectrum shows the lower amplitude A 0 and higher amplitude S 0 
modes which match well with the theoretical dispersion curves. 

In the delaminated plate, the delamination damage, which splits the originally bonded full laminate into two parts 
(above and below the delamination as illustrated in Figure lb), will result in waves propagating in each section. Since 
the two parts, above and below the delamination, will have different layups and thicknesses compared to the original 
pristine plate, waves measured at an outer plate surface region above (or below) the delamination will have modified 
propagation characteristics and wavenumbers. As expected, the delaminated plate spectrum (Figure 5b) shows new f-k 
components in between the original A 0 and S 0 modes, including both positive (forward propagating) and negative 
(backward propagating) counterparts. We compared the new f-k components with the dispersion curves for the material 
above the delamination, equivalent to a two ply [0\0] laminate (the corresponding dispersion curves are plotted in 



Figure6), since only the data of the top plate surface was used for analysis. The comparison in Figure 5b shows, the new 
f-k components match well with the theoretical curves of S 0 mode in the two ply [0\0] laminate. 
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Figure 5. The f-k spectra of EFIT in-plane v y components: (a) for the pristine plate; (b) for the plate with a 
delamination. Solid and dotted lines are theoretical dispersion curves for S 0 and A 0 wave modes, respectively. 




Figure 6. Dispersion curves for a 2 ply laminate with [0/0] layup: (a) f-k curves for 90° direction; (b) group velocity 
curves for 90° direction. 
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Figure 5. The f-k spectra of EFIT out-of-plane v z components: (a) for the pristine plate; (b) for the plate with a 
delamination. Solid and dotted lines are theoretical dispersion curves for S 0 and A 0 wave modes, respectively. 

3.3 Out-of-plane component v z analysis 

Figure 5a plots the pristine plate f-k spectrum for the EFIT out-of-plane v z wavefield. The pristine spectrum shows 
higher amplitude A 0 and lower amplitude S 0 modes which match well with the theoretical dispersion curves. When the 
delamination is introduced, however, new f-k components (positive and negative) appear in the delaminated plate 
spectrum (Figure 5b), and match well with the theoretical curves of A 0 mode in the two ply [0\0] laminate. The analysis 
results of the out-of-plane v z wavefields further confirm that the new f-k components appear in the spectrum of the 
delaminated plate, and are correlated to the waves propagating in the delamination region. 

4. EXPERIMENTAL VERIFICATION 

In sections 2 and 3, the EFIT data show interesting results such as delamination induced wave interactions, 
“trapped waves” and new f-k components. To verify these observations from the EFIT data, this section presents 
comparisons between the out-of-plane EFIT data and experimental data acquired using a ID scanning laser Doppler 
vibrometer (SLDV). 

4.1 Experimental setup 

Two 380 mm x 380 mm 8-ply IM7/8552 composite laminate plates with layup [0 2 /90 2 ] s (the same layup and 
material as used in the EFIT model) were fabricated at NASA Langley. One is a pristine plate and the other is a 
delaminated plate. The delamination is generated by inserting a 20 mm* 20 mm Teflon film under the second ply during 
the ply lay-up process. Figure 6 shows the experimental setup. Guided waves in the plates are excited by a surface 
bonded PZT actuator (APC 851: 7 mm diameter, 0.2 mm thickness) located at the coordinate (0, 0). The excitation 
signal, a 300 kHz 3 -cycle Hanning window smoothed sine wave, is generated by an arbitrary waveform function 
generator (Hewlett Packard 33120A). The signal is amplified to 40V by a voltage amplifier (Krohn-Hite 1506). A SLDV 
(Model: Polytec PSV-400-M2) is used to acquire high spatial resolution velocity wavefield data experimentally. This 
Polytec model is a ID SLDV, and therefore only the motion along the laser beam can be measured. In the test, the laser 
beam is set normal to the surface of the specimen such that the out-of-plane velocity is measured (experimental v z ). 
Scanning is performed along the dotted line (shown in Figure 6) from 5 mm to 60 mm away from the PZT actuator with 
a spatial resolution of 0.2 mm. The SLDV data are used to validate the simulation results in terms of the out-of-plane 
velocity v z . 



Figure 6. The schematic of laser vibrometry experimental test setup. 

4.2 Experimental results 

Figure 7 a and b show the time-space wavefields of SLDV measurements (v z ) for the pristine and delaminated 
plates, respectively. It is seen the A 0 mode has higher amplitude than the S 0 mode. The experimental wavefield of the 
delaminated plate shows waves propagating back and forth (as standing waves) between the boundaries of the delamination 
and exhibit the “trapped waves” phenomenon. The experimental wavefields (Figure 7) verify the observations from EFIT 
out-of-plane wavefields ( 


Figure 4), such as delamination induced wave interactions and “trapped waves” phenomenon. 

Figure 8 a and b show the f-k spectra of the SLDV measurements for the pristine and delaminated plates, 
respectively. Again, the pristine spectrum shows higher amplitude A 0 and lower amplitude S 0 modes. Note that the EFIT 
results in figures 4 and 7 showed even lower amplitude for the S 0 mode. It is expected that this difference between EFIT 
and SLDV data is due to the idealized excitation used in EFIT compared to the real excitation in the SLDV case. When 
the delamination is introduced, new f-k components (positive and negative) appear in the delaminated plate spectrum 
(Figure 8b), and match well with the theoretical curves of A 0 mode in the two ply [0\0] laminate. The analysis results of 
the SLDV measurements are in overall agreement with the results of the EFIT data. The experimental f-k spectra verify 
the new f-k components observed in EFIT f-k spectrum for the delaminated plate. 
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Figure 7. The time-space wavefields measured by the SLDV: (a) for the pristine plate; (b) for the plate with a 
delamination. 
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Figure 8. The f-k spectra of SLDV measurements: (a) for the pristine plate; (b) for the plate with a delamination. Solid 

and dotted lines are theoretical dispersion curves for S 0 and A 0 wave modes, respectively. 


5. DELAMINATION DETECTION AND QUANTIFICATION 

The f-k analysis results of both EFIT simulations and SLDV measurements show new f-k components appear in the 
spectra of the delaminated plate, and are correlated to the waves propagating in the delamination region. This section 
studies delamination detection and quantification by filtering the new f-k components. 

5.1 f-k filtering 

The f-k spectra show that various wave modes are discernible and separable in the f-k domain. Thus, through 
filtering in the f-k domain, any desired wave mode in the f-k spectrum can be extracted [23 , 36, 39]. The f- k filtering 
process can be mathematically expressed as the product between the f-k spectrum U (/, k) and a f-k filter function 
W(J,k) as: 


U w (f,k) = U(f,k)W{f,k) (2) 

where W(f k) denotes the 2D filter in the f-k domain and Uw(f k) is the filtered spectrum. The filtered spectrum Ujff, k) 
is then transformed back into the time-space domain by using the inverse 2D FT in order to obtain the filtered time-space 
wavefield Ujft, x ): 

u w (t,x) = F 2 ^[U w (f,k x )\ (3) 

To extract the waves in the delamination region, the f-k filter W(f k) should be designed as a 2D bandpass filter which 
can keep the desired f-k components correlated to the waves in the delamination region. 

5.2 Delamination detection by using in-plane v y component 

Figure 9a gives the spectral representation of the 2D bandpass f-k filter with a Gaussian profile, which is intended 
to keep the new f-k components, while removing all the rest in the EFIT v y f-k spectrum (Figure 5b). After f-k filtering, 
the bandpass filtered spectrum with the new f-k components is retained, as shown in Figure 9b. The filtered spectrum is 
then inversely transformed back into the time-space domain to reconstruct the wavefield. The reconstructed wavefield is 
presented in Figure 9c, showing that it represents waves propagating forward and backward between the delamination 
boundaries. The filtered wavefield, which only contains the waves in the delamination region, provides an easy means to 
detect and quantify the delamination damage. Details of f-k filtering can be found in reference [39]. 
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Figure 9. Delamination detection by using EFIT in-plane v y component: (a) the f-k filter; (b) the filtered f-k spectrum; 
(c) the filtered wavefield. 


5.3 Delamination detection by using out-of-plane v z component 

Since the out-of-plane data are also available from the SLDV tests, the f-k filtering is applied to both EFIT and 
SLDV out-of-plane data. Figure 10 a, b and c show the f-k filter, filtered spectrum and reconstructed time-space 
wavefield of the EFIT out-of-plane v z component. The reconstructed wavefield shows only the trapped waves 
propagating in the delamination region. Thus, from the reconstructed wavefield, we can detect and quantify the 
delamination location and size. In this case, from the wavefield in Figure 10 c, the delamination is detected aty=30 mm 
with the width 20 mm. The detection result is in consistent with the real delamination location. 

Figure 1 1 a, b and c show the f-k filter, filtered spectrum and reconstructed time-space wavefield of the SLDV data. 
The reconstructed wavefield shows only the trapped waves in the delamination region. From the reconstructed wavefield 
in Figure lie, the delamination is detected aty=30 mm with the width 20 mm. The detection result agrees well with the 



real delamination location and further confirms the feasibility of using f-k filtering for delamination size quantification 
(for this analysis, which is ID in space, the size is quantified in the y-direction length). 
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Figure 10. Delamination detection by using EFIT out-of-plane v z component: (a) the f-k filter; (b) the filtered f-k 
spectrum; (c) the filtered wavefield. 
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Figure 11. Delamination detection by using SLDV measurements: (a) the f-k filter; (b) the filtered f-k spectrum; (c) the 
filtered wavefield. 


6. CONCLUSIONS 

This paper presents our studies of guided wave interactions with delamination damage, as well as delamination 
quantification techniques for laminated composite plates. The 3D EFIT simulations were implemented to acquire both 
in-plane and out-of-plane wave motions for guided wave interactions with delamination damage. Both the in-plane and 
out-of-plane wavefields show waves bouncing back and forth between the delamination boundaries, and distinctively 
exhibit the “trapped waves” phenomenon. The observed behavior agrees with the wavefield observations presented by 
Sohn et al [27]. The wavefields were further analyzed by using f-k analysis to unveil the implicit wave propagation 
characteristics such as f-k components and wave mode contents. In the f-k spectra, new f-k components were observed 
when the guided waves interact with the delamination damage. It is also found that these new f-k components are 
correlated to the waves propagating in the layers above the delamination. 

In addition to the EFIT simulations, the SLDV experimental tests were performed. The SLDV wavefields verify 
the delamination induced wave interaction and “trapped waves” observed in the EFIT out-of-plate wavefields. 
Additionally, the SLDV f-k spectra verify delamination induced new f-k components observed in the EFIT out-of-plane 
f-k spectra. 

A delamination detection method using f-k filtering was studied. By filtering the new f-k components, the waves in 
the delamination region can be extracted. The filtered wavefield provides an easy means to detect and quantify the size 
of delamination damage. The delamination detection method was successfully applied to both the EFIT and SLDV data. 
Though this study was for an idealized simple delamination, it is expected that the f-k analysis can be extended to 
characterize delamination damage in more complex, and realistic scenarios. 
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